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CN ■ Abstract. 

^ ^. The degree of linear polarization in solar flares has not yet been precisely determined 

U ■ despite multiple attempts to measure it with different missions. The high energy range 

^0 ' in particular has very rarely been explored, due to its greater instrumental difficulties. 

QQ , We approached the subject using the Reuven Ramaty High Energy Spectroscopic Imager 

' (RHESSI) satellite to study 6 X-class and f M-class flares in the energy range between f 00 
keV and 350 keV. Using RHESSI as a polarimeter requires the application of strict cuts 

^— j . to the event list in order to extract those photons that are Compton scattered between 

^ ' two detectors. Our measurements show polarization values between 2% and 54%, with 

OO , errors ranging from 10% to 26% in Icr level. In view of the large uncertainties in both 

■ the magnitude and direction of the polarization vector, the results can only reject source 
models with extreme properties. 

0\ ; 

^ , Keywords: Sun; polarization, magnetic field, flares. X-ray polarization 
^ . 

O ■ 

o ■ 

, Measurements of the linear X-ray polarization in solar flares can provide 

' essential information needed to identify the processes responsible for the 

•'-j ■ acceleration of particles and the emission of radiation. For photon energies 

\ in the hard X-ray region, polarization is produced either through electron 

■ Bremsstrahlung or Compton scattering. The degree of polarization is usually 
a complex function of the strength and topology of the magnetic field. Ther- 
mal distributions of electron velocities result in small polarization values 
produced mainly by the conduction-driven anisotropy of the electrons in 
the primary source (Emslie and Brown, 1980) with some contribution from 
photons backscattered in the photosphere. Higher polarization levels are 
expected from non-thermal anisotropic distributions of electrons that are 
accelerated in well ordered magnetic fields. In addition, for such cases the 
observed polarization degree is related to the photon directivity (i.e. the 
anisotropy of the emitted radiation), which depends both on the electron 
beaming details and on the viewing angle. 



1. Introduction 
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Starting in the late seventies, several non-thermal models of X-ray emis- 
sion from solar flares were developed (see McConnell et al., 2002 and ref- 
erences therein). Generally the authors assume a uniform magnetic field 
perpendicular to the solar surface and electrons being accelerated towards 
the chromosphere (Elwert and Haug, 1970; Haug, 1972; Brown, 1972; Bai 
and Ramaty, 1978; Zharkova, Brown, and Syniavskii, 1995). These electrons, 
spiraling downwards along the magnetic field lines, produce Bremsstrahlung 
radiation in collisions with hot plasma. The polarization in the emitted X- 
rays is a function of the energy spectra of the electrons, their pitch angles 
(angle between their velocity vector and the magnetic field), and the column 
density distribution of the ambient plasma. Studying the spectral charac- 
teristics of the detected X-rays provides information about the electron 
energies, while the polarization is a very sensitive tool to sample the electrons 
pitch angle distribution. Both for very small pitch angles, corresponding to 
high electron beaming, and for very large ones, the polarization degree can 
be equally high. The two cases can be distinguished by their polarization 
direction: parallel and perpendicular to the magnetic field line, respectively. 
The predicted polarization values can reach up to 60% at energies above 50 
keV, being even higher at low energies (Haug, 1972). After introducing more 
realistic pitch angle distributions and taking into account photon backscat- 
tering processes in the photosphere, the expected polarization is reduced 
down to 20% or 30% (Brown, 1972; Bai and Ramaty, 1978). 

Several more complex non-thermal models have been developed as well. 
In one of them the magnetic field structure was defined as a semicircular loop 
anchored in the chromosphere (Leach and Petrosian, 1983). This approach 
allowed studying the X-ray emissions separately from different parts of the 
loop. The highest polarization could be produced at the top (up to 85%), 
while the photons observed from the foot-points (in the region of the dense 
chromosphere) would be polarized to the level of around 20%. 

In general, lower energies are predicted to yield stronger polarization 
signals (Haug, 1972; Bai and Ramaty, 1978; Leach and Petrosian, 1983), 
although the strength of this relation can vary depending on the model. In 
some recent theories an opposite trend has also been reported (Zharkova, 
Brown, and Syniavskii, 1995). Experimental verification is usually difficult 
because the low-energy part of the spectrum is strongly contaminated by a 
non-polarized thermal emission. 

In all the cases the observed value of flare polarization is strongly depen- 
dent on the viewing angle. The highest polarization values are expected for 
large angles of view, when the line of sight is perpendicular to the magnetic 
field line. Thus, most theories predict higher polarization for flares located 
near the solar limb. Similar behavior is also expected for the directivity of 
the flare emission: the intensity of the emitted radiation should depend on 
the angle of view (Li et al., 1994; Li, 1995). Constraints related to the model 
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assumptions favor two possible directions of the polarization vector: either 
parallel to the plane defined by the magnetic field lines and the line of sight 
or perpendicular to it (Bai and Ramaty, 1978; Leach and Petrosian, 1983; 
Zharkova, Brown, and Syniavskii, 1995). 

Contrary to the intense theoretical work, only a few polarization mea- 
surements have been conducted in hard X-rays. In this energy range, the 
commonly used technique is based on Compton scattering (Lei, Dean, and 
Hills, 1997). 

First attempts were done at energies around 15 keV by Tindo et al. 
(1970), Tindo et al. (1972), and Tindo, Shuryghin, and Stcffcn (1976) using 
polarimeter instruments on board of several Intercosmos satellites. Although 
initial results (Tindo et al., 1970) showed a linear polarization of around 
40% ± 20%, the later studies of different flares (Tindo et al., 1972; Tindo, 
Shuryghin, and Stcffcn, 1976) found their polarization degrees always com- 
patible with zero. Moreover, the data suffered from limited photon statistics 
and systematic errors related with the detector calibration. 

In measurements of solar flare polarization with the Reuven Ramaty High 
Energy Spectroscopic Imager (RHESSI), Compton scattering can occur in 
a specially installed beryllium scatterer. This method can only be used at 
low energies (20—100 keV). Its details, together with RHESSI's polarimetric 
features, are described in McConnell et al. (2002) and McConnell et al. 
(2004). 

The latest measurements at energies below 100 keV have been performed 
with the SPR-N instrument on board of the Coronas-F satellite by Zhitnik 
et al. (2006). From a sample of 25 solar flares, these authors determined 
upper limits on the polarization degree in the range from 8 to 40% {3a). 
Only for the single case of the 29 October 2003 flare they found a significant 
polarization degree, which increases from about 50% at energies 20 — 40 keV, 
up to more than 70% for the energy channel 60 — 100 keV. 

Recently published, the only two measurements of solar flare polarization 
at high energies (0.2-1 MeV), show la values of 21% ± 10% and -11% ± 5% 
for one flare close to the limb and another near the Sun center, respectively 
(Boggs, Coburn, and Kalemci, 2006). In their studies, the authors applied 
the polarimetric capabilities of the RHESSI satellite in the coincidence mode, 
i.e. without using the Be-scatterer. 

In this paper, we present results in the energy range from 100 keV to 350 
keV, obtained for seven solar flares (X and M classes) also selected from the 
RHESSI instrument database. We used a method based on the scattering 
of photons from detector to detector of RHESSI that has previously been 
applied for polarization studies of gamma ray bursts (Coburn and Boggs, 
2003; Wigger et al., 2004; Rutledge and Fox, 2004). 

In §2 we explain how RHESSI can be used as a Compton polarimeter 
to study linear X-ray polarization at energies > 100 keV. The flares se- 
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lected for analysis, and the criteria used for their selections is detailed in 
§3. Monte Carlo simulations were performed to calculate the response of 
the instrument, and their results are discussed in §4. The final polarization 
results are described and compared with previous measurements in §5. The 
interpretation of the results is done in §6 by comparing with theoretical 
predictions. Finally, a brief summary of the conclusions of our work is given 
in S7. 



2. Method description 

2.1. Instrument 

The Reuven Ramaty High Energy Spectroscopic Imager (Lin et al., 2002) 
was designed to observe solar flares from 3 keV up to 17 MeV. It is capable 
of making spatially, spectrally and temporally resolved images of the Sun 
(Hurford et al., 2002) using the rotation modulation principle (Schnopper, 
Thompson, and Watt, 1968; Skinner and Ponman, 1995). The spacecraft 
is rotating with a period T 4 seconds. The RHESSI angular position in 
the solar coordinate system is calculated using the satellite roll angle, which 
is continuously monitored by the spacecraft aspect systems (Fivian et al., 
2002; Hurford and Curtis, 2002). 

The RHESSI spectrometer (Smith et al., 2002) consists of 9 cooled Ger- 
manium detectors which are split in a thin front and a thick rear segment. 
Low-energy photons arc mostly stopped in the front segments while high- 
energy photons can pass through and reach the rear segments. The energy 
resolution is in the order of a few keV. The arrangement of the 9 detectors 
in the spectrometer is sketched in Figure 1. 

Each photon recorded by RHESSI is characterized by its arrival time, de- 
posited energy, detector number and segment. These parameters are stored 
in the RHESSI event list. The time resolution of RHESSI is equal to one 
binary microsecond (Ib/xs = 2~^° sec). 

Polarization measurements are possible using photons that are Compton- 
scattered from one detector into another one, making a signal in both of 
them. For such events, the effective area is very small, as most of the photons 
are either completely absorbed (in photopeak reactions) or scattered off the 
spectrometer (Hajdas, Wigger, and Zehnder, 2005). Therefore, only a small 
percentage of all photons observed by RHESSI produces more than one 
entry in the event list. The fact that RHESSI rotates allows reducing the 
systematic errors in the polarization analysis. 
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Figure 1. RHESSI detectors seen from the Sun. The grey lines indicate the possible 
scattering directions in the RHESSI fixed coordinate system (§2.4). 

2.2. COMPTON COINCIDENCES 

The Compton scattering probability of photons on free electrons is given by 
the cross section (Klein and Nishina, 1929): 

where tq is the classical electron radius, is the angle between photon infall 
and outfall direction, and r] is the angle between the infall polarization and 
the outfall direction. E and E' are the energies of the photon before and 
after the scattering, respectively. The cross section (u) is minimal if r/ = 
0° and maximal if = 90°. The dependence of the cross section on 77 is 
most pronounced if ^ = 90°. As the RHESSI axis is pointing to the Sun, the 
detectors plane is perpendicular to the photon infall direction, and therefore 
90° for detector-to-detector scattering. 

In order to extract Compton scattered photons from the RHESSI event 
list, several cuts were applied (see Table I). 
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Table I. Selection of cuts introduced to extract Compton scattered photons from 
the event list. 



Criteria 



Values 



Time interval 
Coincidence width 
Single event energy 
Coincidence energy (sum) 
Kinematical cut 
Detector segments 
Pairs selection 
Multiplicity 



Between 1 min. and 4 min. around the flare peak 
1 b/us 

25 keV <Ei< 300 keV 

100 kcV <Ei + Ej < 350 keV 

6* = 90° ± 45° 

Rear segments (without detector 2) 
Neighbor detectors 
Two-event coincidences only 



As most of the high energy photons pass through the front segments 
without any interaction, only the rear segments were used for the polariza- 
tion analysis. This condition ensures that Compton scattering events happen 
within the same b//s (Wigger et al. (2004), Equation (8)). Detector 2 was not 
taken into account as it operates in a different way than the others (Smith 
et al., 2002). 

When more than two events happen at the same time, information about 
the polarization direction is lost. Therefore we selected photons which were 
registered in exactly two detectors, calling them 'coincidences'. Furthermore, 
only coincidences between neighboring detectors were chosen, due to the very 
low probability of direct scattering between remote detectors (see Figure 1). 

The sum of the two energies of the coincidence {Ei and Ej) had to fall 
into our selected range between 100 and 350 kcV. In addition, the energy cut 
imposed on the individual detectors was in the range 25-300 keV. Below this 
region, the detector noise and the level of accidental coincidences (see §2.3) 
increased very strongly. A final constraint was applied by introducing the 
so-called kinematical cut which excludes coincidences that are incompatible 
with Compton scattering, using the kinematical relation between the scat- 
tering angle 6 and the observed energies. The kinematical cut was especially 
useful to reject photons scattered from the Earth atmosphere (see §4). 

Applying the above cuts, we obtain a raw list of coincidences which is 
still contaminated by background events of different origin. 

2.3. Background subtraction 

Two major sources of background were taken into account in the follow- 
ing analysis. The first are accidental coincidences which occur when two 
independent solar photons are simultaneously detected. The second kind of 
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background is not related to the flare itself but produced by cosmic rays and 
the cosmic gamma ray background. 

The rate of accidental coincidences is proportional to the square of the in- 
coming photon flux. In order to obtain the number of accidental coincidences 
(-^acc) we repeated the same procedure that was used to flnd the raw list of 
coincidences (Table I), but taking pairs of events which are time delayed by 
20 b/xs < dt < 30 hjis (see Wigger et al., 2004). As the typical detector dead 
time is in the order of several microseconds, the delay length was chosen to 
be longer than this. The rate of accidental coincidences over this range is 
approximately constant and was used for backgroTind subtraction. 

The second kind of background was determined either before or after the 
flare. For this purpose we selected time intervals shifted by 18 orbits, i.e. 24 
hours, before or after the flare peak. In this way, the spacecraft geomagnetic 
coordinates were similar (sec first part of Table II for an example) , and the 
systematic effects coming from the background variations along the orbit 
were strongly reduced. Proper subtraction of this background requires some 
additional conditions: no other flare or high energy event should be present in 
that period, and the spacecraft operational status regarding attenuator state 
and decimation logic should be close to the one in the moment of the flare 
observation. The total number of background coincidences obtained in this 
way is denoted by Nfyg^to^. The period selected for background subtraction 
contains a (usually negligible) number Ni,g^acc of accidental coincidences. It 
has to be subtracted from Nf^g^t^t because it is already included in Nacc cal- 
culated during the flare peak. The number of coincidences produced by the 
non-flare related background is then Ni,g = Ni^g^tot ~ ^bg,acc- As an example, 
the time evolution of the background signal observed one day before the 
20 January 2005 flare is displayed together with the flare in Figure 2. It 
reproduces very well the background levels observed before and after the 
flare peak. 

Table II. Geomagnetic coordinates of RHESSI satellite and numbers of coincidences determined 
at the peak of the flare on 20 January 2005 and during the background measurement period. 
The data collecting time was 240 seconds. 



PEAK OF FLARE 



BACKGROUND 



Starting time 



06:43:00 20-Jan-2005 06:42:59 19-Jan-2005 



RA(°) 
Dec(°) 



-76.07 



37.91 



-80.26 



37.91 



Total coincidences 



43313 ± 208 



5874 ± 77 



Accidental coincidences 
Solar flare Compton scattering coincidences 



26907 ± 35 



10598 ± 225 



66 ±2 
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Figure 2. Signal and background coincidences vs. time for the flare on 20 January 2005 
together with the background levels measured approximately one day before at the same 
position. 



Finally, the rate of Compton Scattering events was calculated for any 
time t according to 

Nc{t) = Ntot{t) - Naccit) - Nkg{t + r), (2) 

where T is chosen such that both the geomagnetic coordinates and the 
orientation (roll angle) of RHESSI are in the background period as close as 
possible to the ones during the flare peak. N^ot is the number of coincidences 
in the flare peak according to Table I, Nf,g is the corresponding number of 
background events from flare-unrelated photons, and Nacc is the number of 
accidental coincidences. 



2.4. Polarization analysis 



For each coincidence, the line connecting the two detectors involved in the 
Compton scattering process defines the angle used to construct the modu- 
lation curve. According to the RHESSI fixed coordinate system, there are 
four possible scattering directions d (labeled 0°,45°,90°, and 135°) between 
neighboring detectors (Figure 1). In each direction, the number N^iji) of 
coincidences per time bin (with ti = U MOD T) was divided by their sum 
in the same direction (nd(Ti) = v-^tw'* J - Such normalization is needed in 

order to eliminate the dependence on the number of detector pairs which is 
different for each direction (see Figure 1). 

The RHESSI roll period T was measured using the spacecraft roll aspect 
system (Fivian et al., 2002; Hurford and Curtis, 2002) which provides the 
RHESSI angular position with an accuracy of 1 arcmin. We obtained the 
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satellite rotation period for each flare using a linear fit. In order to verify its 
stability and check for possible drifts, T was calculated in 1 sec steps over a 
720 sec interval containing the flare. The maximum variations between the 
measurements were below 1 ms, thus negligible comparing with the bin size 
used in the polarization analysis (r^ 167 ms). It was also found that the drift 
of the period is smaller than (2.1 ± 2.0) • 10~^. It corresponds to a change in 
the rotation period by less than 0.05 ms, giving an upper limit of 0.53° (Scr) 
on the possible phase shift for the longest time interval used in the analysis. 
The asymmetry curve can be constructed as follows: 



^0-90 (Tj) = -— -— . (3) 

no(rj) + ngoiTi) 

The coincidences occurred in the directions 45° and 135° were properly 
shifted and included in the equation above. Using asymmetry to deter- 
mine the polarization minimizes the effect of the lightcurve variations and 
grouping to have only two directions improves the statistics. 

In order to relate the time variable r with angular directions in the Sun, 
we used the relation: T]i = —^-Ti + rjo. By convention, 77 = corresponds 
to solar West, and rj = it/2 corresponds to the solar North, r/o is the an- 
gular position of RHESSI X-axis (0° direction from Figure 1) with respect 
to the solar West in the moment when measurement started. Trough this 
coordinate transformation, we obtained ^o-9o(''") in heliocentric coordinates: 

Due to the sinusoidal dependence of the Compton cross section for the 
scattering of a polarized photon (Equation (1)), the asymmetry curve is also 
a sinusoidal function with period equal to a half of the RHESSI rotation. It 
can be represented by a function: 



^(?7) =/Xp-cos(2(r?-<^ + 7r/2)), (4) 

where the amplitude jip is a positively deflned value equal to the flare 
modulation factor. The phase 4) is the polarization angle from the flare. 

Comparing the experimental amplitude jjLp with the modulation factor 
/Ltioo from Monte Carlo simulations for a 100% polarized flux (§4), allows to 
determine the polarization degree of the solar flare (H) (McConnell et al., 
2002): 

H = (5) 
Afioo 

Two parameters, H and ^, fully describe the polarization state of a solar 
flare and are needed for comparison with theoretical predictions. 
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Figure 3. Distribution of the flares in the solar disc. Their position was extracted from 
images taken with RHESSI. 



3. Flares selection 



The following criteria were applied to select the flares for polarization anal- 
ysis: large intensity, strong high energy component and negligible contam- 
ination with particles (cither from the flare itself or from the radiation 
belts). Also, since theory predicts highest polarization close to the solar 
limb (Brown, 1972; Bai and Ramaty, 1978; Leach and Petrosian, 1983), we 
focused on limb-close flares. After applying all these conditions, we were able 
to select six X and one M class flares, five of them located within less than 
120 arcsec from the solar limb (Figure 3). 

In order to study polarization in the most explosive part of the energy 
release, only the peak of the flare was chosen (see Figure 4). The time period 
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Lightcurve of Flore on 23-Jul-2002 
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Figure 4- Ligthcurve of the solar flare on 23 July 2002 and background selected for it. 

The number of single events registered in rear segments of all detectors (except number 
2), with energies between 25 keV and 300 keV is plotted as a function of time. The part 
of the flare used for the polarization measurement is enclosed by the two vertical dashed 
lines. 



for analysis varied between one and four minutes depending on the duration 
of the flare peak. 

The spectrum of each flare was analyzed to find the energy ranges of 
different emission mechanisms. For this purpose, a fit was performed, using 
the RHESSI OSPEX fitting tool (Tolbert, 2006), with a combination of a 
thermal Bremsstrahlung curve and a broken power law (Figure 5). In all 
cases, the thermal contribution was found to be negligible at energies above 
50 keV. For further analysis we selected the non-thermal Bremsstrahlung 
region (100-350 keV). The spectral indices are given in Table III and corre- 
spond to a single power law fit of that part of the flare. These values were 
subsequently used in the simulations performed to determine the instru- 
mental response function and its polarization modulation factor (see §4). 
Although photons with energies between 50 and 100 keV are already in the 
non-thermal emission region, their interaction in the RHESSI detectors is 
governed by photoelectric absorption, leaving only a marginal number of 
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Fitted Spectrum of Flore on 19-Jon-2005 
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Figure 5. Background-subtracted spectrum of the flare on 19 January 2005. The fit was 

made for energies between 12 keV and 600 keV combining a thermal curve and a broken 
power law. Thermal component is negligible above 40 keV, and the background becomes 
dominant around 500 keV. 



Compton scattering events. On the high-energy side, the threshold value 
was chosen to avoid regions dominated by background. 

The main parameters that describe the flares selected are summarized in 
Table III. 
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Flare number (RHESSI) 


2072301 


3110221 


4111002 


5011710 


5011911 


5012005 


5082502 


Date 


23 Jul. 2002 


2 Nov. 2003 


10 Nov. 2004 


17 Jan. 2005 


19 Jan. 2005 


20 Jan. 2005 


25 Aug. 2005 


Start time 


00:18:00 


17:03:00 


01:59:00 


09:35:36 


07:57:20 


06:21:24 


04:33:48 


End time 


01:15:44 


18:00:36 


02:36:52 


10:38:48 


09:03:32 


07:27:04 


04:55:56 


Duration (s) 


3468 


3456 


2272 


3792 


3972 


3940 


1328 


Class 


X4.8 


X8.4 


X2.6 


X4.0 


X1.4 


X7.1 


M7.0 
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60 


150 


240 


240 


Counts (s"^) (E > 25 keV) 


67421 


126440 


28296 


50001 


20798 


128520 


21804 


Spectral index (100-350 keV) 


3.1 


3.5 


3.4 


3.8 


2.9 


3.2 


3.2 
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Figure 6. Simulated asymmetry plot for a 100% polarized flare with spectral index -3.1. 

4. Simulations 

Monte Carlo simulations have been performed to calculate the response of 
the RHESSI polarimeter to a 100% polarized solar flare. For this purpose, 
the exact mass model of the whole satellite and its germanium spectrometer 
has been constructed and implemented in the GEANT 3.21 simulation 
code (CERN, 1994). The spacecraft was illuminated by a uniform beam 
of photons coming parallel to the RHESSI rotation axis. Each incoming 
photon was fully tracked and all energy depositions made on its path in 
any of RHESSI detectors were recorded. In this way a simulated event list 
was created and subsequently used to determine the modulation factors. 
The procedure to extract the modulation factors from the event list was the 
same as for the analysis of the solar flares (§2), but neglecting the accidental 
coincidences. We made two sets of simulations: 

1. The modulation factors for 100% polarized emissions were determined 

for each flare using a photon energy distribution with spectral indices 
as given in Table HI. The incoming photons were 100% polarized and 
their energy range was either 100-600 keV or 80-350 keV. It was found 
that the contribution from photons above 350 keV to the modulation 
curve in our energy range (100-350 keV) was negligible. An example of 
a modulation curve is shown in Figure 6. 

2. In the second case we took monoenergetic photons with 100% polar- 
ization at thirteen energies between 100 keV and 1000 keV (Figure 7). 
The mean modulation factor of each flare was subsequently calculated 
by averaging the monoenergetic modulation factors. The weights used 
for the average were proportional to the number of coincidences per 
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Exponential fit of 100% modulation factor 
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Figure 7. Simulated energy dependence of jUioo- The dashed line corresponds to the fit 
function: Atioo = 56 • e-^/''«° - 8668 • e-^/^^ 



energy bin. /xioo is maximum around 170 keV. It decreases for lower 
energies due to the low energy threshold (25 keV) which only allows 
detecting recoil electrons from photons scattered at very large angles. 
At such angles, modulation factor of the instrument is very small, in 
accordance to Equation (1). Above 170 keV hiqq diminishes with energy 
following the polarization sensitivity based on Equation (1) for photons 
scattered around 90°. Using two exponential functions with properly 
chosen coefHcients reproduces the above features in a simple way (see fit 
function in the capture of Figure 7). 

Both approaches gave the same results. For example, for the case of the 
23 July 2002 flare, with an spectral index equal to -3.1, we obtained with the 
first method /xioo = 32.8±1.6% and with the second one /xioo = 32.4ib5.4%. 

The analysis of the solar flare lightcurves from single detectors revealed 
periodic structures that can be attributed to photons reaching RHESSI after 
being scattered by the Earth's atmosphere. The contamination of the mod- 
ulation curves by such photons was computed with another set of Monte 
Carlo simulations. For this purpose a simplified system consisting of the 
satellite and the Earth with its atmosphere was constructed. The Earth was 
represented by a solid sphere with twelve layers of atmosphere extending 
up to about 50 km above the surface. The mass as well as the chemical 
composition of all the atmospheric sheets were equal, while the density varied 
in accordance with their height. 

Simulations of Earth scattering were performed using unpolarized pho- 
tons with spectral indices typical for the analyzed flares and their corre- 
sponding angular positions between RHESSI, the Earth and the Sun. The 
largest fraction of photons detected from the atmosphere was found when the 
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Figure 8. Simulations of scattering by the Earth's atmosphere. The figure shows single 
events with energies between 100 keV and 350 keV arriving at the rear part of detector 
number 4. The solid line shows the total photon flux, while the flux coming directly from 
the Sun appears as a dashed line. 



Earth reached the angle of 90° with respect to the RHESSI-Sun direction. 
In the energy range 100 keV-350 keV, up to 30% of the observed photons 
were coming from the Earth, producing a strong modulation in the single- 
event lightcurves (Figure 8, solid line). The influence of such photons on the 
asymmetry curves, extracted from coincidences, was much smaller. Most of 
the Earth scattering caused only accidental coincidences or did not pass the 
kinematical cut. Finally, the contamination by the Earth-scattered photons 
was, in the worst case, less than 8%. Considering their low level and flat 
distribution along the modulation curves, the modifications of the measured 
modulation factors were negligible compared to the overall statistical error. 

The presence of the grids above the detectors could produce periodic 
structures also in the coincidence lightcurves and therefore this effect was 
carefully studied with Monte Carlo simulations. Those sources which are 
not in the center of the Sun are, during RHESSI rotation, intermittently 
obscured by the grids. This produces in their single-event lightcurves a 
modulated profile with 180° periodicity (Hurford et al., 2002) and such a 
pattern could be mixed up with the real polarization signal. Firstly, the 
lightcurve modulation is strongly reduced for high energy photons reaching 
the rear detectors. Secondly, the width of the angular bin used for polariza- 
tion analysis is large enough to average all count-rate variations caused by 
the grid modulation. The lightcurve variations further cancel out since the 
grids of detector pairs used for coincidences are generally aligned in different 
directions. 
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We performed simulations of a flare situated in the solar limb where the 
grid effect would be the strongest. Fine grids have been approximated keep- 
ing the slit/sIat width ratio while the coarse ones were exactly implemented. 
The modulation factor obtained for a non-polarized flux of 10*^ photons in 
the 100-350 keV energy range and with a power law spectrum index of 
-3.1, was 2.5% ± 1.9%. This is similar to the signal measured from a non- 
polarized source situated in the center of the sun, where grids do not cause 
any modulation. Therefore, the grids effect can be neglected. 

Notice that due to the positive definiteness of the polarization degree, 
even an unpolarized signal gives a non- vanishing amplitTidc in the asymme- 
try plot when applying to our analysis. From a 0% polarized simulated flare 
with around 8500 coincidences, we obtained a modulation factor equal to 
(3.4 ± 1.6)%. 

5. Results 

5.1. RHESSI POLARIZATION MEASUREMENTS 

The asymmetry curves A{r]i) are displayed in Figure 9, together with the 
best fit of the function in Equation (4) . These curves have a periodicity of 
T/2 (T equal to RHESSI rotation period). To improve the statistics, the 
second half of the asymmetry curves was added to their first half, plotting 
only the range 0°-180°. 

The degrees of polarization H were calculated using Equation (5) with 
/xioo taken from simulations. The mean value of /xioo was 33.2% and its vari- 
ations between different flares were below 1.1%. The resulting polarization 
degrees were found to be between 2% and 54%, with error bars varying from 
10% to 26% at the la level (see Table IV). The H values of all the flares are 
plotted versus the flare class in Figure 10, where no significant correlation 
between these two quantities can be observed. 
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Figure 9. Asymmetry curves of all flares analyzed, extracted for photon energies between 
100 keV and 350 kcV. The thick lines show the best fits with the function from Equation 
(4). The angle of the minimum in the fit curve indicates the flare polarization direction 
in heliocentric coordinates. 



solflarepol3.tex; 5/02/2008; 13:26; p. 18 



Table IV. Summary of polarization results for the flares studied, /ip is the observed modulation factor. 11 is the polarization degree of 
the flare, and (j> its polarization angle given in heliocentric coordinates. 



Flare number (RHESSI) 


207230f 


3110221 


4111002 


5011710 


5011911 


5012005 


5082502 


Date 


23 Jul. 2002 


2 Nov. 2003 


10 Nov. 2004 


17 Jan. 2005 


19 Jan. 2005 


20 Jan. 2005 


25 Aug. 2005 


Ntot 


7439 ± 86 


34723 ± 186 


3816 ± 62 


2142 ± 46 


5688 ± 75 


43313 ± 208 


6139 ± 78 


Naca 


2269 ± 10 


21427 ± 31 


506 ±5 


473 ±5 


783 ±6 


26907 ± 35 


602 ±5 


Nbg 


1758 ± 42 


5135 ± 72 


2047 ± 45 


733 ± 27 


2784 ± 53 


5808 ± 77 


3717 ±61 


Nc 


3411 ± 97 


8160 ± 202 


1262 ± 77 


937 ± 54 


2121 ± 92 


10598 ± 225 


1820 ± 100 


Atloo (%) 


33.0 ±1.6 


32.4 ± 1.8 


32.9 ± 1.7 


33.4 ± 1.9 


31.6 ± 1.5 


34.4 ± 1.6 


34.4 ± 1.6 


Hp (%) 


0.6 ± 4.5 


9.1 ±3.9 


11.7 ±8.6 


9.4 ±8.2 


17.1 ±6.5 


7.2 ±3.3 


2.2 ± 8.4 


<P (deg) 


151 ± 195 


96 ±12 


104 ± 24 


71 ±29 


170 ± 11 


66 ±14 


102 ± 104 


n m 


2 ± 1 1 
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Figure 10. Results on the degree of polarization, with their la error bars. 



The data does not show any preferential direction of polarization. When 
plotting the polarization angles of the flare sample in heliocentric coordinates 
(Figure 11, left), the points seem to concentrate around the North-South 
solar direction but the error bars are too large to extract a firm conclusion. 
When the angle between the polarization direction and the line that joins 
the flare position and the center of the Sun was calculated, no tendency was 
found (Figure 11, right). 

Despite the very high single-photon count rates of more than 20000 counts 
per second, the mean number of Compton events found per flare was only 
around 4000, and the smallest one does not even reach 1000 coincidences. 
This reflects the small value of the RHESSI effective area for Compton 
polarimetry as discussed in §2. The signal-to-background ratio is on the 
average only around 0.5. Depending on the flare, the largest background 
contribution is produced either by accidental coincidences or by cosmic 7- 
ray background. The number of counts found with each flare are compiled 
in the first part of Table IV, where the errors refer to la and are purely 
statistical. 
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Figure 11. Polarization degree vs. polarization angle plotted in two diflFerent reference 
systems: with respect to solar equator {left) and with respect to the radial line that joins 
the center of the Sun and the flare position {right). 



5.2. Comparison with other measurements 

With respect to the polarization amplitude, our results are consistent with 
previous measurements made at higher energies by Boggs, Coburn, and 
Kalemci (2006). In the particular case of the 23 July 2002 flare that they 
also analyzed, our value (2% it 14%) is smaller than theirs (21% it 10%), 
but agreement is found at the 1.5(7 level. The difference can be explained by 
the different time periods and energy ranges used in both cases. Extending 
the time and energy windows for our analysis towards the values selected by 
Boggs, Coburn, and Kalemci (2006) provides very similar polarization levels. 
The polarization angles of the two flares measured by Boggs, Coburn, and 
Kalemci (2006) are aligned along the North-South direction in heliocentric 
coordinate system. However, these authors conclude that polarization is az- 
imuthal for near-limb flares, but radial for those close to the Sun center. We 
can not confirm this rule from our observations (see Figure 11, right). The 
values found in the present work are more uniformly distributed between 35° 
and 85°, independently of the flare location. Taking into account the values 
of the error bars and the size of the statistical sample of analyzed flares, 
further measurements aimed to verify the observed disparity are needed. 

Recent polarization data at energies up to 100 keV have became available 
from measurements with the SPR-N instrument on board of the Coronas-F 
satellite (Zhitnik et al., 2006). From a sample of 25 solar flares, upper limits 
on the polarization degree were found to be in the range from 8 to 40% 
(3o"). These values are in good agreement with our results (typically within 
2a level). In particular, for the flare on 20 January 2005, observed simulta- 
neously by both satellites, the polarization value from RHESSI observations 
was equal to 21%± 10% while the upper limit in the Coronas-F measurement 
was equal to 17%. Again, a more direct comparison is not possible because 
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both the energy range and the time intervals analyzed were different. For one 
flare (on 29 October 2003), the Coronas-F instrument showed a significant 
polarization degree that increases from about 50% at energies 20-40 kcV, 
up to more than 70% for the energy channel 60-100 keV. Unfortunately, 
RHESSI polarization analysis of the 29 October 2003 flare was not possible 
due to a high contamination of its detectors with charged particles. 



6. Interpretation 

Our results were compared with theoretical predictions for the non-thermal 
photon emission given by Bai and Ramaty (1978), Leach and Petrosian 
(1983) and also with the 0% polarization hypothesis. Bai and Ramaty (1978) 
provide the most comprehensive set of theoretical data, covering the whole 
energy range from 10 to 500 keV. These authors considered primary X-ray 
emission due to Bremsstrahlung of the accelerated electrons moving towards 
the photosphere, adding also the X-ray Compton backscattering component. 
The photon polarization was studied for different electron spectra and pitch 
angle distributions and results were presented as a function of the observing 
angle. At energies around 200 keV, typical for our analysis, the polarization 
reaches maximum values between 20 and 30% and decreases quickly for 
flares located closer to the solar center. Taking this trend into account, the 
predicted polarization degrees for our sample of flares range from -19% to 
1%, where the negative sign indicates a polarization direction parallel to the 
magnetic field, and the positive one perpendicular to it. 

Leach and Petrosian (1983) analyzed emissions from more complex, loop- 
shaped magnetic fields, but most of their polarization prediction is given only 
for two energies: 16 and 102 keV. As the dependence of the polarization on 
the energy was weak, we extrapolated their values towards our energy range. 
The highest polarizations for large energies are expected to come either from 
the top of the fiare loop, or from the transition region above the chromo- 
sphere. The latter option roughly corresponds to RHESSI observations in 
which the high-energy emissions come from the foot-points. Nevertheless, 
a clear distinction between chromospheric and transition zone emission was 
not possible. Leach and Petrosian (1983) proposed several models depending 
on the magnetic field gradient, the electron pitch angle distribution and the 
spectral index. We selected three different cases for comparison purposes: 
one with a homogeneous magnetic field and a uniform distribution of the 
electron pitch angles (model 3) , another with equally homogeneous magnetic 
field but pitch angles close to 90° (model 5), and the last one with a large 
magnetic field gradient and large pitch angles (model 8). After correcting 
for the flare position on the Sun, the expected polarization degrees of our 
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Table V. Reduced from all the flares combined, 

obtained by comparison of several theories with our 

measurements. 

MODEL Reduced 

Bai and Ramaty (1978) 1.03 

Leach and Petrosian (1983) (model 3) 0.86 

Leach and Petrosian (1983) (model 5) 2.70 

Leach and Petrosian (1983) (model 8) 0.83 
0% polarization 0.82 



flares are in the ranges between -7% and 18% for model 3, from -45% to 
-80% for model 5, and from -6% to 23% for model 8. 

The latest theoretical work on polarized emission, with the most advanced 
electron beam dynamics, was presented by Zharkova, Brown, and Syniavskii 
(1995). Polarization calculations were done only up to 100 kcV (for photon 
spectral indices typical to our flares) , and showed increased values for higher 
energies. Their predictions, extrapolated to the energy range used in the 
present work, are similar to those from Bai and Ramaty (1978) and cannot 
be distinguished by the following analysis. 

The comparison of our observations with the models was done by gen- 
erating asymmetry curves for the flares polarized in accordance with the 
theoretical predictions. Monte Carlo simulations were performed for each 
flare individually, using exactly the same RHESSI mass model as in §4. 
The expected polarization angle and amplitude were properly determined 
by adjusting the theoretical predictions to the flare position on the Sun. The 
high number of simulated events allowed to keep the statistical errors on a 
very low level comparing with the experimental data. Finally, the reduced 

value was calculated using the data points from all seven flares as given 
by Equation (6). 

1 n m (K[e _j^t\2 

^ n{m - 1) {rif^^ {^Nl,f + {^Nl^f ' ^ ^ 

where i is the flare number and j is the angular bin in the asymmetry plot. 
A^Pj and Nlj are the experimental and theoretical numbers of coincidences 
respectively, and AiV?^- and AA|,^- are their statistical uncertainties. The 
total number of degrees of freedom was equal to 77. The results for all four 
models and the 0% polarization hypothesis are displayed in Table V. 

Prom the values of Table V we can reject the model 5 from Leach and 
Petrosian (1983), with a 90% of confidence. For the rest of the models the 
X^ are very close to unity, preventing us from distinguishing between them. 
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Within the error bars, they all agree with the experimental data equally 
well. The same is also valid for the 0% polarization hypothesis. Further 
refinement would require a polarimeter able to deliver data with error bars 
on the level of 1-2 percent. 



7. Summary and conclusions 

Measurements of hard X-ray polarization have been performed for six X- 
class and one M-class flares from the RHESSI database. Our flare sample was 
identified after applying selection criteria to the signal strength, background 
levels and level of contamination by charged particles in RHESSI detectors. 
The selected energy range, from 100 to 350 keV, connects the old (Tindo 
et al., 1970; Tindo et al., 1972; Tindo, Shuryghin, and Steffen, 1976) and 
new (Zhitnik et al., 2006) results at low energies with the only measurement 
reported at a high energy band (Boggs, Coburn, and Kalemci, 2006). 

We found values for the polarization degree in the range between 2% and 
54%, with statistical errors from 10% to 26% at the la level. The polarization 
angles are distributed between 66° and 170° in heliocentric coordinates. 
They do not show any preferential orientation of the polarization, neither 
parallel nor perpendicular, with respect to the radial direction defined by 
the position of the flare in the Sun. In addition, no significant dependency 
between the orientation of the polarization and the distance of the flares to 
the Sun center was found. 

The polarization orientation with respect to the line that joins the two 
major foot-points of the flare was also studied. For this purpose, images of 
the flares were constructed with RHESSI at different energies. The emis- 
sion above --^30 keV was found to be produced around the foot-points. 
However, no correlation between polarization direction and foot-points ori- 
entation was found. Similarly, no relationship between solar flare intensity 
and polarization degree could be observed. 

The polarization degree from the 23 July 2002 flare measured by Boggs, 
Coburn, and Kalemci (2006) at high energies is in agreement with our results 
at the 1.5(7 level. Their conclusion about the orientation of polarization re- 
spect to the radial direction passing trough the flare position (perpendicular 
to it for flares in the limb and parallel for flares close to center) cannot 
be confirmed by our measurements. Our angles are distributed between 35° 
and 85° independently on the flare location. Comparison with theoretical 
predictions is more complex, as the direction of the polarization is expected 
by the theory to change around 300 keV (Bai and Ramaty, 1978). Again, 
more observations with better accuracy are needed. 

Regarding the low energy measurements, our data are in good agreement 
with the recent results from the SPR-N instruments on board of the Coronas- 
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F satellite (Zhitnik et al., 2006), typically within 2(j. Their large sample of 
25 solar flares reveals low polarization degrees, providing 3a upper limits 
from 8 to 40%. Unfortunately, no information is given about the polarization 
direction. In the 29 October 2003 flare, for which the Coronas-F team claims 
a polarization of 70%, RHESSI measurements were contaminated by charged 
particles. 

Theoretical predictions of the non-thermal Bremsstrahlung emission pro- 
vide polarization levels of the order of 20%. However, depending on the 
assumptions used by different authors, the expected polarization can differ 
not only by its value, but also by its orientation (see, for example, the 
calculations at 100 keV from Bai and Ramaty (1978) and Leach and Pet- 
rosian (1983)). As the statistical uncertainties provided by our instrument 
are in the same order as the model predictions, equally good agreement is 
found for any of them as well as for the case of 0% polarization. In order to 
distinguish between different models, polarimetry measurements of at least 
2% accuracy are needed. Only model number 5 from Leach and Petrosian 
(1983), predicting polarization up to 85%, could be rejected by our RHESSI 
data. In this model, the magnetic field has the same strength at the top and 
the bottom of the loop, and the electrons spiral at pitch angles close to 90°. 

RHESSI has made the first steps towards the understanding of the po- 
larization phenomena in solar flares above 100 keV, where non-thermal 
emission dominates. Measurements with accuracy better than 10 to 20% 
were, however, hardly possible. This is due to its small effective area and 
high levels of flare-induced background. Continuation of such studies will 
require a dedicated polarimeter that must solve the problems inherent to 
the RHESSI design. Emphasis should be put on increasing the effective area 
and improving the background rejection capabilities. In particular, a better 
time resolution will reduce the number of accidental coincidences, and the 
optimization of the detector dimensions will improve the detection efficiency 
for the Compton scattering. 
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